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The gauge theory of electroweak interactions contains 
a striking feature. In quantum electrodynamics, the pho- 
tons carry no electric charge and thus lack photon-to- 
photon couplings and do not self- interact. In contrast, 



the weak vector bosons carry weak charge and do in- 
teract amongst themselves through trilinear and quartic 
gauge boson vertices. 

The most general WW-y and WWZ interactions can 



be described |l|, l2| using a Lorentz invariant effective La- 
grangian that contains fourteen dimensionless couplings, 
seven each for WW-f and WWZ. Assuming electromag- 
netic gauge invariance and CP conservation reduces the 
number of independent couplings to five (electromagnetic 
gauge invariance requires gj — 1), and the Lagrangian 
takes the form: 



wwv = igYiW^WV - W£V V W'* V ) 
gwwv 

+ iKvWZw v V>" + %*wlw» v v* x 

where denotes the W boson field, W llv =df 1 W v — 
d v Wn, = dfj,V v - dvVjj,, F=7 or Z, and M w is the 
mass of the W boson. The global coupling parameters 
gwwv are gww-y—— e and gwwz = —ecot6w> as in the 
standard model (SM) in which e and 9w are the magni- 
tude of the electron charge and the weak mixing angle, re- 
spectively. In the SM A 7 =A^=0 and g1=gf =ft 7 = Kz=l. 
For convenience, anomalous trilincar gauge couplings 
(anomalous TGCs) Any and Agf are defined as Ky — 1 
and gf — 1, respectively. 

The W boson magnetic dipole [iw and electric 
quadrupole qw moments may be expressed in terms of 
the coupling parameters as 

g 



m w 

As mentioned above, g^=^- 

If the coupling parameters have non-SM values then 
the amplitudes for gauge boson pair production grow 
with energy, eventually violating tree-level unitarity. The 
unitarity violation can be controlled by parametrizing the 
anomalous couplings as dipole form factors with a cut- 
off scale, A. The anomalous couplings then take a form 
a(s) = ao/(l + s/A 2 ) 2 in which y/1 is the center-of-mass 
energy of the colliding partons and ag is the coupling 
value in the limit •§ — [3j . The quantity A is physically 
interpreted as the mass scale where the new phenomenon 
responsible for the anomalous couplings is directly ob- 
servable. The cutoff A is conservatively set at the limit 
of sensitivity, close to the collision center-of-mass energy. 
We use A = 2 TeV; coupling limits depend only weakly 
on A for A > 1 TeV in hadronic collisions at Tevatron 
energies. 

We measure the electroweak coupling parameters 
through the study of gauge boson pairs. Several pro- 
cesses contribute to SM boson pair production. Fig.QJa) 
shows i-channel production of dibosons in which V1V2 are 
WW, WZ, or W7. The s-channel production shown in 
Fig. DJb) involves boson self-interactions through a tri- 
linear gauge vertex. Final states (V1V2) produced via the 
WW Z coupling are WW or WZ. Final states produced 



4 




FIG. 1: Vector boson pair production via (a) f-channel and 
(b) s-channel diagrams. For Vi = W and V2 = "t/Z, Vo = W. 
For Fi = V 2 = W, V = j/Z. 



through the WWj coupling are WW or VF7. The typical 
effect of anomalous TGCs is to increase the cross section 
especially at high boson transverse momentum (pt)- We 
thus analyze corresponding observables to measure such 
effects. 

Previously published limits on anomalous TGCs from 
a combination of channels come from the DO Collabora- 
tion in the 1992-1996 Tevatron run with integrated lumi- 
nosity (£) of 100 pb" 1 [4], the CDF Collaboration with 
the current Tevatron run (£-350 pb- 1 ) [5], and LEP2 
experiments [6] . The best previously published W boson 
magnetic dipole moment result is from a combination of 
measurements by the DELPHI Collaboration [7j. 

In this Letter, we investigate the IFVF7 and WWZ tri- 
linear vertices through diboson production. We set limits 
on the non-SM or anomalous TGC parameters Ay, Any, 
and Agf . These limits are derived from a combination 
of previously published measurements involving four final 
states: W^iv-f @], WW->£v?v 0, WZ^lvl't' [l^, 
and WW jW Z^ivjj [ll|, in which i is an electron or 
muon, v is a neutrino, and j is a jet. Each measure- 
ment used data collected by the DO detector [l2j from 
pp collisions at y/s = 1.96 TeV delivered by the Fermilab 
Tevatron Collider. 

The process IF7— >£i/j is sensitive only to the WW^f 
coupling. The process was studied with data corre- 
sponding to 0.7 fb _1 0. The main requirements were 
an electron with transverse energy Et>25 GeV or a 
muon with transverse momentum pt>20 GeV, a pho- 
ton with Et>9 GeV, missing transverse energy ^ T >25 
(20) GeV for the electron (muon) channel, and separa- 
tion between the photon and lepton in 77— 4> fl3| space 
of AR=^J(Arf) 2 + (Ac/)) 2 >0.7. Furthermore, to sup- 
press final state radiation the three-body transverse mass 
|14| of the lepton, photon, and $ T was required to ex- 
ceed 120 (110) GeV for the electron (muon) channel. 
In total 180 (83) candidate evj (^7) events were ob- 
served. After subtracting backgrounds, the signal was 
130± 14 sta t±3.4 sys t (57±8.8±1.8) events, consistent with 
the SM prediction of 120±12 (77±9.4) events for the ev^ 
(fivj) channel. The photon Et spectra of the Wj can- 
didates in the data and those estimated for the back- 
grounds are input into the combination. For W7 pro- 
duction in the presence of TGCs, spectra were simulated 
using the Baur Monte Carlo (MC) [H, Ell with a fine 
grid in A 7 — Ak 7 space. 
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The WW—*lvt'v measurement [9j used data corre- 
sponding to an integrated luminosity of 1 fb _1 . The data 
were divided into three channels defined by the flavor of 
the leptons from the W boson decays: ee, efi, and 
For all channels, the leading lepton had pt>25 GeV and 
the trailing lepton had pt>15 GeV. The leptons were re- 
quired to have opposite charge. In the data 22 (ee), 64 
(e/i) and 14 (/U/i) candidate events were observed, con- 
sistent with the sum of SM WW and backgrounds of 
23.5±1.9 (ee), 68.6±3.9 (e/x) and 10.8±0.6 (m) events. 
Two-dimensional histograms of leading and trailing lep- 
ton pt were produced for the data and backgrounds and 
used as inputs in the combination. Distributions for SM 
and anomalous TGC values were generated using the 
WW/WZ event generator from Hagiwara, Zeppenfeld, 
and Woodside (HZW) @. 

The WZ^lvt'i measurement pjj selected the four 
final states eee, ee/x, /i/ie, and The data corre- 

sponded to an integrated luminosity of 1 fb _1 . All three 
charged leptons were required to have p-r>15 GeV. Z bo- 
son candidates consisted of like-flavor lepton pairs with 
mass 71<M ee <lll GeV or 50<M^<130 GeV. For the 
eee and fififi channels, the oppositely charged lepton pair 
with mass closest to the Z pole mass was chosen as the 
Z boson candidate. To select W boson candidates, the 
fi T must have exceeded 20 GeV. To reduce background 
events from ti to a negligible level, the magnitude of the 
vector sum of the charged lepton transverse momenta and 
the $ T was required to be less than 50 GeV. The sum 
over all channels yielded 13 candidate events in the data 
consistent with a SM estimate of 9.2 ± 1.0 W Z events 
and 4.5 ± 0.6 background events. The pt of the Z bo- 
son is sensitive to anomalous TGCs and is used in the 
combination. The HZW MC is used to estimate the SM 
spectrum as well as spectra from anomalous TGCs. 

Finally, the WW/W Z—^lvjj measurement [ll[ se- 
lected events in which one W boson decays leptonically 
and the other boson decays hadronically. The data 
corresponded to an integrated luminosity of 1.1 fb _1 . 
The main requirements were an electron or muon with 
Pt>20 GeV, $ T >20 GeV, and at least two jets with 
Pt>20 GeV with the leading jet satisfying pt>30 GeV. 
In total 12,473 (14,392) candidate events in the evjj 
(fJ-vjj) channel were observed, consistent with the SM 
prediction of 12,460±550 (14,370±620) evjj (pvjj) 
events [l7| ■ An observable sensitive to anomalous TGCs 
is the pt of the dijet system. The data and background 
spectra for this variable are used as inputs for the com- 
bination. Spectra with anomalous TGCs were generated 
with the HZW MC. 

Distributions of the sensitive observables mentioned 
above for each final state are generated for signal with 
the corresponding Monte Carlos and for backgrounds us- 
ing simulations or data. The signal distributions vary 
as a function of the TGC parameters under study both 
in spectral shape and event yield. In addition to allow- 
ing variation in the TGC parameters themselves, nui- 
sance parameters are used to allow systematic offsets to 



vary within their uncertainties. A simultaneous fit to the 
data distributions is performed in order to determine the 
anomalous TGC limits. The \ 2 function used in this fit 
is El: 



X = 



ffi C p {d i]mi {R)) C G (R k a k -Q,a k ) 
~ 2in HI £^;^ 11 



fc=i 



£ G (0;0,a k ) 



= 2 mi (R) — di — di In 



rrii(R) 



k=l 



in which the variables i and k index the number of his- 
togram bins (Nb) and the number of systematic uncer- 
tainties (N s ) respectively. In this function £ p (a;(3) is 
the Poisson probability for a events with a mean of j3 
events; C G (x; h,<j) is the Gaussian probability for the 
value x in a distribution with a mean value of \i and a 
variance a 2 ; R k (in vector form as R) is a dimensionlcss 
parameter describing departures in nuisance parameters 
in units of the associated systematic uncertainty a k ; di 
is the number of data events in bin i; and rrii(R) is the 
number of predicted events in bin i. The number of bins 
used in the fit is the sum of the number of bins in each 
kinematic distribution for each channel. 

In total 49 sources of systematic uncertainty are con- 
sidered. As implied in Eq.[T] systematic uncertainties are 
treated as Gaussian priors on the expected number(s) of 
events. Systematic uncertainties on the luminosity, lep- 
ton identification, and theoretical uncertainties on the 
cross sections for the backgrounds estimated from MC 
are correlated across all observables. Uncertainties on 
background estimates based on data are correlated across 
specific final states within a diboson production channel 
as appropriate. The uncertainties with the largest im- 
pact on the result are those related to background cross 
sections and the luminosity. The effect of incorporating 
systematic uncertainties into the fit is to degrade the re- 
sulting limits by ~30%. 

Four two-dimensional surfaces in TGC space are exam- 
ined: (a) each of the three pairings of the three free pa- 
rameters (Ak 7 , A, Agf ) while respecting SU (2) l®U 
symmetry by using the constraints Akz = A<?f — 
AK 7 tan 2 9w and Xz — A 7 = A [n| and (b) the (Ak, A) 
plane for the equal-couplings scenario in which k 7 = 
kz = k, A 7 = Xz = X. The two-dimensional 68% and 
95% C.L. contours are shown in Figs. and [3] The 
two-dimensional contours for W boson magnetic dipole 
and electric quadrupole moments are shown in Fig. 3] 
The one-dimensional 68% and 95% C.L. limits for each 
coupling parameter, with the other couplings parameters 
fixed at their SM values, are shown in Table HI 

These results provide the most stringent limits on 
anomalous values of WWj and WWZ TGCs measured 
from hadronic collisions to date. The 95% C.L. lim- 
its in both scenarios represent an improvement relative 
to the previous DO [H and CDF Q results of about a 
factor of 3. When respecting SU(2)l <8> U(l)y symme- 
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TABLE I: One-dimensional \ 2 minimum and 68% and 95% 
C.L. allowed intervals on anomalous values of WW*y and 
WWZ TGCs. Note that (jlw and qw are in units of (e/2Mw) 
and (e/M w ) respectively. 



Results respecting 


SU{2) L ®U{1) Y 


symmetry 


Parameter 


Minimum 


68% C.L. 


95% C.L. 


Ak-, 


0.07 


[-0.13,0.23] 


r n oo n qqI 
[— u.zy, u.ooj 


Ao? 


0.05 


[-0.01,0.11] 


[-0.07,0.16] 


A 


0.00 


[-0.04,0.05] 


[-0.08,0.08] 


AW 


2.02 


[1.93,2.10] 


[1.86,2.16] 


qw 


-1.00 


[-1.09,-0.91] 


[-1.16,-0.84] 




Results for equal-coupling 


s 


Parameter 


Minimum 


68% C.L. 


95% C.L. 


Ak 


0.03 


[-0.04,0.11] 


[-0.11,0.18] 


A 


0.00 


[-0.05,0.05] 


[-0.08,0.08] 


fJ.W 


2.02 


[1.94, 2.09] 


[1.88,2.15] 


qw 


-1.02 


[-1.09,-0.94] 


[-1.16,-0.87] 



try, our measurements with 68% C.L. allowed intervals 
of k 7 - 1.07±°;|g, A = 0.00±°;°i and gf = 1.0518$ 
are only factors of approximately 2-3 times less sen- 
sitive than the combined results from the four LEP2 
experiments: k 7 = 0.973l£j^, A = -0.028l£jg? and 
gf = 0.984toml, also at 68% C.L. <o\. Furthermore, 
with only 1 fb _1 of data our sensitivity is comparable to 
that of an individual LEP2 experiment @, S [HI, HI . 

We also extract measurements of the W boson mag- 
netic dipole and electric quadrupole moments. When 
respecting SU(2)l <8> U(l)y symmetry with gf=l we 



measure 68% C.L. intervals (one-dimensional with the 
other parameter held at its SM value) of = 
2m±° f 9 (e/2M w ) and q w = -1.00 ± 0.09 (e/M w ), 
respectively. The most stringent previously pub- 
lished result is fj>w — 2.22l ; 19 (e/2Mw) and qw = 
-1.18lg;^ (e/M w ) from the DELPHI Collaboration 0). 

In summary, we presented measurements of anoma- 
lous WW j and WWZ trilinear gauge couplings and re- 
lated W boson magnetic dipole and electric quadrupole 
moments based on the combination of four diboson pro- 
duction and decay channels using 0.7—1.1 fb _1 of data 
collected with the DO detector at the Fermilab Tevatron 
Collider. While many of the measurements considered 
in this combination are limited by statistics, projections 
indicate that a combination of CDF and DO data with 
5 fb _1 each will improve the sensitivity to levels compa- 
rable or better than the combined LEP2 limits. 
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FIG. 2: Two-dimensional 68% and 95% C.L. limits when re- 
specting SU (2) l ® U(1)y symmetry and assuming A = 2 TeV, 
for (a) A vs. Ak 7 , (b) A<?f vs. Ak 7 , and (c) A<?f vs. A. In 
each case, the third coupling is set to its SM value. 
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FIG. 3: Two-dimensional 68% and 95% C.L. limits for A vs. 
Ak when enforcing the equal-couplings constraints and as- 
suming A = 2 TeV. 




FIG. 4: Two-dimensional 68% and 95% C.L. limits for the 
W boson electric quadrupole moment vs. the magnetic dipole 
moment (a) when respecting SU(2)l <8> U(1)y symmetry and 
(b) when enforcing equal-couplings constraints. In both cases 
we assume A = 2 TeV. 



